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The Lazarus project was designed to make the most of limited 3D binary black-hole simulations, 
through the identification of perturbations at late times, and subsequent evolution of the Weyl 
scalar $4 via the Teukolsky formulation. Here we report on new developments, employing the 
concept of the "quasi-Kinnersley" (transverse) frame, valid in the full nonlinear regime, to analyze 
late-time numerical space-times that should differ only slightly from Kerr. This allows us to extract 
the essential information about the background Kerr solution, and through this, to identify the 
radiation present. We explicitly test this procedure with full-numerical evolutions of Bowen-York 
data for single spinning black holes, head-on and orbiting black holes near the ISCO regime. These 
techniques can be compared with previous Lazarus results, providing a measure of the numerical- 
tetrad errors intrinsic to the method, and giving as a by-product a more robust wave extraction 
method for numerical relativity. 
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I. INTRODUCTION 

The strong-field interaction of black-hole binary sys- 
tems — from early approach through capture, mutual 
orbit and eventual merger, to ring-down of the end-state 
single hole — is expected to be a primary source of 
gravitational radiation at all frequency scales, and has 
been a focus of theoretical and numerical attention for 
forty years. Early perturbative studies and two- 

dimensional numerical evolutions of axisymmetric bina- 
ries (head-on collisions) •3, '3, '3, '3, were successful 
in producing late-stage waveforms representing gravita- 
tional radiation. However, the move to full 3D simu- 
lations of more general initial-data sets has proved ex- 
tremely difficult. Evolutions using the "standard ADM" 
3-1-1 decomposition of Einstein's equations and simple 
zero-shift gauge conditions have stable lifetimes of ~ 
15 — 30M (where M is the total mass of the space-time), 
far too short a time to complete a useful physical simula- 
tions, much less extract the gravitational radiation emit- 
ted — the ultimate aim of numerical source simulations. 

The Lazarus project 0,0] was conceived in the con- 
text of such limitations. Working under the assumption 
that the late-evolution 3-f 1 data can be considered a per- 
turbation of a single Kerr black hole, Lazarus extracts the 
radiation content everywhere in the numerical domain, 
and uses it as initial data for a Teukolsky perturbative 
evolution. In this manner, the original simulation may 
be extended almost indefinitely, long enough to capture 
the entire development of the outgoing radiation. 

The beginning of the full- numerical simulation can also 
be interfaced with a Jar limit approximation method and 
similar techniques to evaluate a common regime of ap- 
plicability can be developed j^. Here, for the sake of 
definiteness, we will assume a set of initial data as provid- 
ing this interface values and focus on the full-numerical 
/ close limit matching. 



Lazarus has been very successfu l, p roducing the first 
convergent wave forms IT^ IT^ [l5l | from 3D evo- 
lutions. In most "plateau" was identified — a 
range of extraction times T where the emitted energy 
remained flat and consistent. This plateau begins when 
the 3-1-1 data is linearly perturbed from Kerr, and should 
end only when the radiation has begun to leave the 3-f 1 
numerical domain altogether; at this latter time, Teukol- 
sky extraction will no longer capture the full radiation 
content. However, the numerical instability of the 3-1-1 
simulation may pollute the Teukolsky initial data for late 
extraction times. Somewhat surprisingly, such a plateau 
seems to exist even in situations where a common appar- 
ent horizon has yet to form, e.g., short-lived evolutions of 
so-called "ISCO" [Innermost Stable Circular Orbit) and 
"pre-ISCO" runs. 

Meanwhile, great strides have been made in full 3D 
simulations over the last decade, due to the casting of 
Einstein's e qua tions into more numerically stable formu- 
lations [IliillillilllSilllllllillillilli^ 

|2^ . the development of advanced techniques for han- 
dlin g th e singularities inherent in black-hole space-times 
[23. I29I l30| . and the availability of increased compu- 
tational resources, coupled with mesh-refinement tech- 
niques [H m m in III El . The culmination of these 
advances is several successful evolutions of black-hole- 
binary systems past the symbolic "one orbit" barrier 
[37Ll3a |. For physical systems requiring less than ~ lOOM 
of evolution time to reach a quiescent final state, these 
improvements enable the direct extraction of radiation 
from the 3-1-1 fields, whether through Weyl curvature 
components or Regge- Wheeler- Zerilli variables 39, 40]. 

These advances, however, do not lessen the relevance 
of the Lazarus project. Despite great progress in recent 
years, it is fair to say that the problem of black-hole 
binaries is not completely solved. The most astrophysi- 
cally interesting simulations can still not be evolved long 
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enough to reach their assumed quiescent state. Directly 
extracted radiation is still calculated at observer loca- 
tions that lie in the "near-field zone" , or may be under- 
resolved at more distant locations, and is in general pol- 
luted by poor outer-boundary conditions. As long as such 
limitations exist, there is a place for perturbative meth- 
ods such as Lazarus. Besides, numerical simulations are 
still very computationally intensive, and avoiding the last 
Ri lOOM of binary black-hole evolutions means saving 
days to weeks of supercomputer time. 

However, Lazarus makes approximations in its ap- 
proach. Principal among these is the set of ad hoc choices 
needed to translate the 3-1-1 curvature information into a 
Kerr background -I- perturbations. The validity of these 
choices will depend on the data being evolved, and is 
difficult to quantify a priori. 

In this paper, we update the Lazarus project in light of 
recent work on transverse frames, in a way that may help 
resolve some of these issues. Beetle et al. have pro- 
posed a way of identifying the principal directions of a nu- 
merical space-time in the 3-1-1 ADM split. This method 
— local in nature — allows us to narrow the gap between 
the numerical tetrad and the Kinnersley tetrad appropri- 
ate to the Teukolsky evolution without any background 
assumptions. When calculated with such a tetrad, the 
longitudinal Weyl scalars and vE's will vanish, while 
the "monopole" scalar ^2 will take on its Kinnersley- 
tetrad value, and the transverse scalars ^'o and ^'4, which 
carry the radiative degrees of freedom, will differ from 
their Kinnersley-tetrad values only by a complex fac- 
tor. This remaining factor can be compensated for, on 
a known Kerr background in Boyer-Lindquist (BL) co- 
ordinates, via a single spin-boost transformation at each 
point in space. 

Although the quasi-Kinnersley frame by no means re- 
moves all uncertainties from the problem of radiative ex- 
traction, it goes sufficiently far that we expect it to im- 
prove the Lazarus procedure considerably. In particular, 
we expect that the different — and more rigorous — 
path to Teukolsky waveforms will give us an error esti- 
mate for the tetrad dependence of original results, while 
the robustness of the new technique should allow us to 
attempt consistent wave extraction from earlier in a nu- 
merical evolution. 

Additionally, the quasi-Kinnersley frame may achieve 
much in the simpler problem of direct radiation extrac- 
tion ,41^ 1,42, uMj • In the past , we used an ap- 
proximate tetrad to calculate the Weyl scalar. Since the 
quasi-Kinnersley tetrad can be constructed locally, with- 
out knowing Kerr parameters and BL coordinates, we 
should now achieve a better approximation to the Kin- 
nersley tetrad during the 3-1-1 evolution, and thus directly 
extract waveforms without needing the background data. 

The remainder of this paper is laid out as follows: 
in Section m we summarize the essentials of the origi- 
nal Lazarus procedure for constructing Cauchy data for 
Teukolsky evolution, as well as some of the main re- 
sults from this procedure. In Section IlIII we review the 



concepts of transverse frames, and the quasi-Kinnersley 
frame, and describe the numerical implementation of 
these concepts in our evolution code. In Section Hvl we 
present comparative results from the application of old 
and new techniques for three test problems, two of which 
have already been addressed with original Lazarus [isj . 
Discussion of the results obtained, and future work can 
be found in Section Appendix ^ contains expres- 
sions for several quantities related to the evaluation of 
Weyl scalars in a numerical tetrad in BL coordinates; 
Appendix IbI contains perturbative results for the three 
test cases used in Section Hvl 

A. Notation and Conventions 

In the rest of the paper, we shall generally assume 
a metric signature of (— , -|-, Our use of alge- 
braic quantities related to the Kerr-BL solution is non- 
standard, but consistent with pj| : we use an additional 
quantity A ~ r^ + for compactness. 

Our sign convention for the definition of the Weyl 
scalars is such that for the Kerr solution in BL coor- 
dinates (and using the Kinnersley tetrad), the only non- 
vanishing scalar is ^2 — M/{r — ia cos 9)^. 

Vector quantities are denoted by an arrow overhead, 
except for unit vectors, which are instead capped by a 
circumflex Q. Complex conjugation is denoted by an 
overbar (~). 

II. THE ORIGINAL LAZARUS METHOD 

Here we provide a summary of the "original" Lazarus 
procedure; full details can be found in |^ , which we shall 
refer to as Paper I from now on. 

We start with a full- numerical 3-1-1 evolution of a 
space-time of interest. At some extraction (coordinate) 
time T — before the simulation crashes due to numerical 
instabilities — we map the evolved data to a black-hole 
perturbative evolution code p6L |47| . This simpler code 
can then be evolved stably for as long as is needed to 
determine the full history of the gravitational radiation 
generated. 

To assess the level of deviation from Kerr at late times 
in the 3-1-1 evolution, introduced an invariant quan- 
tity, the speciality index: 

S^27J^/I^, (2.1) 

where the two complex curvature invariants I and 
are essentially the square and cube of the self-dual part, 
Cabcd = Cabcd + (V2) Cafcrnn (^""cci' °^ ^hc Weyl tcusor: 

I = Cabcd C^"''" and J = Cabcd C^"^^ C'""''^ (2.2) 

The geometrical significance of S is that it measures 
deviations from algebraic speciality (in the Petrov clas- 
sification of the Weyl tensor) . For the unperturbed alge- 
braically special (Petrov type D) Kerr solution, 5 = 1. 
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However, for interesting space-times involving nontrivial 
dynamics, like distorted black holes, which are in gen- 
eral not algebraically special (Petrov type I), we expect 
5 = 1-1- AS, and the size of the deviation AS ^ 0, with 
leading second perturbative order, can be used to assess 
the applicability of black-hole perturbation theory. 

As the expected end-state of most interesting black- 
hole simulations is a single spinning (Kerr) hole, the per- 
turbative code implements the Teukolsky equation 
The Kerr metric in BL coordinates takes the form: 
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where ip = \E'4 is the spin-2 Teukolsky function, and 
^ = r -\- ia cos 9. Here p = m°' la-b ttt-^ is a Newman- 
Penrose spin coefficient, and is a Newman-Penrose 
Weyl scalar, both calculated using the Kinncrsley tetrad. 

In the Newman-Penrose formalism 50] , there are actu- 
ally five complex Weyl scalars, formed from contractions 
of a null tetrad {1°", n", m", m") with the Weyl tensor: 



(2.5) 



The encode all the vacuum curvature information 
of the Weyl tensor. As space-time scalars, they are 
coordinate-independent; however they do depend on the 
particular null tetrad used. With an appropriate tetrad, 
in weak- field regions, the interpretation of the ^'i is as fol- 
lows: embodies the "monopole" non-radiative gravi- 
tational field; ^'i and '^3 contain the longitudinal radia- 
tive degrees of freedom (ingoing and outgoing, respec- 
tively), while and ^'4 contain the physical transverse 
radiative degrees of freedom (ingoing and outgoing, re- 
spectively) |51|. For a numerical space-time that con- 
tains a Kerr hole plus perturbative gravitational waves, 
^4 should contain only the appropriate outgoing radia- 
tion. 



The asymptotic behavior of solutions to the Teukol- 
sky equation is best expressed in terms of the so-called 
tortoise coordinate r*: 



In 



2M 



In 



2M 



r± = M ± x/M^-a^. 



(2.6) 



The point r» = roughly corresponds to the location of 
the maximum of the Kerr solution's scattering potential 
barrier (see, for example, Eqn (415) and preceding mate- 
rial in Chap. 8 of ^5^). For this reason, it should not be 
crucial to obtain initial data for the Teukolsky equation 
all the way down to the horizon (r = r+ => = —00), 
as long as we have data for some < 0. 

Thus to perform the Teukolsky evolution of radiative 
data that corresponds to the late-time evolution of our 
3-1-1 initial data, we must identify the parameters (M, a) 
of the Kerr background, and calculate the radiative Weyl 
scalar ^4 using the Kinnersley tetrad. Estimation of 
the physical parameters can be performed fairly reliably 
through identification of physical invariants such as the 
apparent horizon and ADM mass or correcting (itera- 
tively) the initial data parameters by the radiative losses. 
Evaluating ^4 using the correct tetrad is less straightfor- 
ward. 



A. Tetrad Choice 

In BL coordinates, the Kinnersley null tetrad takes the 
form ,53.1 : 



iKin = ^[A,A,0,a], 
"Kin = ^ [A, -A,0,a] 



(2.7) 



1 



iriKu 



V2C 



i a sin9, 0, 1, ■ 



smf 



Using this tetrad, the spin coefficient p takes the form: 
p = (r- cos6l)~i = l/C (2.8) 

However, the BL coordinates will not, in general, coin- 
cide with the numerical coordinates used in the full 3-1-1 
evolution. We can address this issue in a post-processing 
step after the evolution, but we must still extract enough 
curvature information to construct Rather than 

output all the components of the Weyl tensor Cabcd, it 
is more efficient to calculate the Weyl scalars with a nu- 
merically convenient tetrad, and transform the results to 
the Kinnersley values during post-processing. 

The simpler tetrad we use during evolution is a sym- 
metric null tetrad constructed from the unit hypersur- 
face normal f and a set of three orthonormal unit spatial 
vectors e(i) = eg, 6(2) = e^, 6(3) = e^, suitably orthonor- 
malized via a Gram-Schmidt procedure: 

^m.ni = ^(7" + 6(3)) , Wnum = ;^(t - 6(3)), 
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"iimm = -^(6(1) + ie(2)). (2.9) 

Similar tetrads have been commonly used in radiation 
extraction from 3+1 numerical investigations 113,0, 
IH^, and such a tetrad was used in the earliest investiga- 
tions of the asymptotic radiative degrees of freedom of 
the Weyl tensor |5J • If we have long-lived 3D numerical 
evolutions, whose physical domain extends far from the 
strong-field region, the '^4 extracted should yield a good 
measure of the actual outgoing gravitational radiation. 
We will refer to (|2.9|l hereafter as the numerical tetrad; 
explicit formulas for the Kerr-BL Kinnersley tetrad are 
given in (|A1|I . 



consist of a combination of type-I and type-II tetrad ro- 
tations parametrized by A, followed by a type-Ill "spin- 
boost" transformation parametrized by the real scaling 
factor Fa and pure-phase complex factor Fb ■ This trans- 
formation carries over to the Weyl scalars; in particular, 

mf"' = [(i:)-l)2^'f;"''^ + 4i^(D-l)*r'" 
-6 -4iA{D + l) ^-^""^ 

+ (i? + l)2^r"] /(4J^1^^|), (2.15) 

where D = \J + 1. The equivalent transformations for 
the other Kinnersley-tetrad scalars can be found in Eqns 

(iX4i-ir7i. 



B. Reconstructing Boyer-Lindquist Coordinates 



D. Summary of the Original Lazarus Procedure 



The reconstruction of the BL coordinates (i, r, (p) is 
highly non-trivial. We approach the problem via the fol- 
lowing ad hoc steps: 

(i) . assume no polar coordinate distortion; 

(ii) . assume that with maximal slicing, numerical time 

approaches Boyer-Lindquist time; 

(iii) . derive the radial coordinate from the equatorial 

value of numerical X; 

(iv) . add a radius-dependent correction to the numer- 

ical azimuthal coordinate that zeros out the off- 
diagonal three- metric component 7^0 . 



In short: 



TBL = {3MVI] 



1/6 



(2.10) 
(2.11) 

(2.12) 



(eBL=7r/2) 

0BL - 0num = 0offsct== / {ir- 4> / lcl>4,) dr' . (2.13) 

oo 

These (last two) coordinate transformations can only 
be performed as a post-processing step, after the termi- 
nation of the full 3D numerical evolution. 



C. Transforming to the Kinnersley Tetrad 

For the Kerr-BL metric, we can move from the nu- 
merical tetrad to the Kinnersley tetrad by a set of linear 
transformations governed by the parameters 



A = a shi9\ 



, Fa 




Fh 



c 



(2.14) 



(Note that \Fb\ = 1, since CC = r^ + a'^ cos^ = S.) The 
specific transformations can be found in Eqn IIA3II : they 



To summarize, the original Lazarus procedure involves 
the following steps: 

(i) . construct at every point on the numerical grid a 

coordinate null tetrad of the form l|2.9|l : 

(ii) . calculate the corresponding Weyl scalars \I>™™, the 

tetrad invariants X, J^, and the speciality index S, 
which we monitor during the entire full-numerical 
evolution; 

(iii) . determine the BL coordinates from the numerical 

ones via the transformations (|2.1(JI - [?.13|I : 

(iv) . use the transformation H2.15|l to obtain the 

Kinnersley-tetrad ^'4 and 9t^'4; 

(v) . evolve the Cauchy data ip = p^^ ^4 and 9tV' using 

the Tcukolsky equation (|2.4|) : 

(vi) . extract the gravitational radiation information, 

such as waveforms and total energy radiated at dif- 
ferent extraction times T. 



E. Results of Original Lazarus 

The original Lazarus procedure as outlined above has 
been applied extensively to numerical evolutions of black- 
hole binary data, extending from head-on collisions, graz- 
ing collisions and putative circular-orbit data at various 
orbital separations 0, 0, 0| . 

One of the main shortcomings of Lazarus is the ad hoc 
nature of the coordinate transformations in H2.10l - [TTS|l . 
Even with data in the linearized regime, Lazarus will be 
sensitive to how well the numerical space-time satisfies 
these coordinate assumptions. For instance, we have no 
guarantee that there is no angular distortion between the 
numerical and BL radial coordinates (although studies 
of post-merger apparent horizons have noted a definite 
tendency of the horizon shape to circularize in numerical 
coordinates). 
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An even bolder assumption is that maximal slicing 
will yield a late-time lapse with the Kerr-BL shape. 
The Kerr-BL lapse satisfies the maximal slicing equation; 
however, the numerical lapse also depends on the bound- 
ary conditions. In practice, we use Dirichlet boundary 
conditions, with values equal to that of the Schwarzschild 
lapse (M = 1) at the same coordinate position. This 
choice will give us a lapse shape (and hence numerical 
time) qualitatively like Kerr-BL; however we only have 
experimental quantitative experience with the quality of 
the fit. It can be shown for instance that the numerically 
obtained lapse quickly takes the form of the BL lapse 
over points in the exterior of the horizon; the greatest 
deviations occur within the "potential barrier" at r*. A 
plot demonstrating this for evolved QCO binary data is 
shown in Fig. 6 of Paper I. 

A complete and unambiguous solution of the coordi- 
nate problem is not yet available. Until it is, it seems 
sensible to work to minimize Lazarus's coordinate de- 
pendence. One obvious area to address is the transforma- 
tion H2.15|l of the curvature from numerical to Kinnersley 
tetrad; the numerical Weyl scalars are both mixed and 
scaled by coordinate-dependent factors. An improvement 
has been made possible by evaluating the Weyl scalars 
instead with the quasi- Kinnersley tetrad, which is the 
subject of the next section. 



III. LAZARUS WITH THE 
QUASI-KINNERSLEY FRAME 

In this section, we review some recent results on trans- 
verse frames in numerical space-times, and discuss how 
they can be used to develop new techniques to improve 
the Lazarus procedure described in Section^ The main 
aim is to move closer to the calculation of the actual Kin- 
nersley tetrad during the full 3-1-1 evolution, and thus to 
minimize the need for ad hoc coordinate and tetrad cor- 
rection schemes. 



A. Transverse Frames and the Quasi-Kinnersley 
Frame 

The Kinnersley tetrad on a Kerr background is trans- 
verse - when calculated with this tetrad, the "longitudi- 
nal" Weyl scalars and vanish. This continues to 
hold for the non-trivial situation of a Kerr hole plus per- 
turbing radiation. However, the coordinate and slicing 
ambiguity of a generic numerical space-time, even when 
only perturbatively different from Kerr, mean that this 
tetrad can be difhcult to identify. The Lazarus tetrad- 
transformation procedure outlined in the previous section 
will not, in general, yield a transverse tetrad. 

A way of identifying transverse tetrads is through find- 
ing the eigen-bivectors of the self-dual Weyl tensor Cabcd 
|41| (see also Chapter 4 of When expressed in 

a 3-1-1 decomposition, the Weyl tensor is projected to 



Cac = Cabcd t''' , and the eigen-equation becomes 

C^a" = {E\-iB\)a'' ^ Xa", (3.1) 

where E"'f^ and B"'^^ are the so-called electric and magnetic 
parts of the Weyl tensor [s^. fSTf . This can be recast by 
projecting the Weyl tensor onto the orthonormal triad 
{e(i)}: following 56] , we can write (summation is implied 
over parenthetical indices): 

Cac = Q(i){j)^(i)a^(j)c^ (3-2) 

a" = y(,)e^,), (3.3) 

where Q(i)(j) is a symmetric complex 3x3 matrix whose 
components are: 

Q(i)(i) = -*2 + (I'd + 1'4)/2, 

Q(l)(2) = 1(^-0 -^'4)/2 , Q(l)(3) =*3-^'l, 
Q(2)(2) = -*2-(*0 + *4)/2, (3.4) 
Q(2)(3) = -i(1'l + *3) , Q(3)(3) = 2^-2, 

and the ^'^ are calculated using the symmetric tetrad 
(|2.9|) . Then the eigen-equation H3.1|) reduces to 

The three eigenvalues A of the matrix Q{i){j) label three 
transverse frames; the corresponding eigenvectors give 
the principal directions of the space-time. One of the A 
is preferred — it will be analytic near the point 5 = 1 
in the complex iS-plane. The frame related to this A 
is called the quasi- Kinnersley frame. The preferred A 
will be numerically equal to 2 x the value of the 

monopole Weyl scalar as calculated with the Kinnersley 
tetrad. 

Beetle et al. ||41| have described the analytic determi- 
nation of the Weyl eigenvalues. In practice, establishing 
analyticity of the the eigenvalues is not necessary numer- 
ically - Mars '5^ has pointed out that close to Kerr, the 
eigenvalue with the largest complex norm will give the de- 
sired frame; this conclusion has been made more secure 
by |4H . who have shown that this is a valid conclusion 
everywhere in the disc ||5 — 1|| < 1. Instead of following 
the analytic route, therefore, we use the LAPACK rou- 
tine zgeev [59I I to determine numerically the eigenvalues 
and eigenvectors of Q{i){j) ■ We select the largest-modulus 
eigenvalue as the appropriate (quasi-Kinnersley) one. 



B. Tetrad Reconstruction 

Beetle et al. 0| lay out a procedure for constructing 
a tetrad (?°, n"', m°, rh"')qK in the quasi-Kinnersley frame 
from the eigenvector V(i): first construct cr" from (|3.3|l : 
normalize it so that || a \\— 1; separate into real and 
imaginary parts: cr" = -t- iy"; construct a third or- 
thogonal vector = s°'ij^x^y'^; then the new null tetrad 
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vectors are: 

HqK = 
TOqK = 



where c = |c|e*-^ is an arbitrary spin-boost parameter. 
We take |c| = l,x = 7''/2, in order to produce a tetrad 
that asymptotes to the original numerical tetrad at large 
distances. The subscript "qK" will refer to this specific 
choice of c from this point. 

If it happens that the eig envector g ° is identically real, 
II X \\~ 1, II y ||~|| z ||~ ^11 X ||2 — 1 = (see Appendix 
Section lBSl for analytic examples in the "close-slow" limit 
of binary Bowen-York data) . This is not a problem when 
constructing /qK and UqK; however Eqn H3.8|l will now 
be undefined. A valid complex null vector mqK can still 
be formed in this case, by replacing the second term in 
parentheses by any linear combination a-\-ih oi two unit 
spatial vectors orthogonal to x. Beetle et al. supply one 
such choice in Eqn (29) of However, the resulting 
TOqK at this point will have real and imaginary compo- 
nents that may not match continuously to neighboring 
points. 

In a continuous domain, we could imagine evaluating 
TOqK by taking a limit from neighboring points; on a nu- 
merical domain (especially in 3D), this is an impractical 
approach, as (i) it would necessitate knowing in advance 
which points would need to be interpolated, and (ii) we 
would need a very dense numerical mesh to carry out such 
an interpolation. Additionally, there are cases where the 
pathological points are not isolated, but cover the entire 
domain. This is the case for Schwarzschild and Brill- 
Lindquist data. In such cases, no interpolation procedure 
is possible. 

For these reasons, we use an alternative tetrad recon- 
struction procedure, one that avoids pathologies entirely. 
We start by following the prescription of Eqns H3.6l - rr7jl 
for the reconstruction of the real null vectors ZqK and 
riqK, as in |4l| : these vary smoothly from point to point, 
even for pathological regions when || z \\ = 0. 

Next, we take the original numerical complex null vec- 
tor and split it into real and imaginary parts: m = 
X -\- lY . Now starting from these vectors, orthonormal- 
ize them according to a set of Gram-Schmidt-like steps. 
Since ZqK and nqK are already correctly orthonormalized, 
the remaining requirements are: 



TOqK • ^qK 
TOqK • ^qK 
TOqK • TOqK 
TOqK ■ TOqK 



O^X-?qK = 0, y.^qK^O, 

^ X • nqK = 0, f • nqK = 0, 
O^X-f = 0, X-X = Y -Y, 
X ■ X + Y - Y ^1. 



The last two equations combine to imply that: 

X ■ X = Y - Y = 1/2. 

To impose these conditions, we begin with X = 
6(1) /\/2, and enforce the conditions in turn (note that 

since /qK and riqK are null, the Gram-Schmidt procedure 
looks slightly unusual): 



X 
X 



X + {X ■ /qK)nqK, 
X + {X ■ nqK)/qK, 



X -> x/y2{x-x). 

In a similar manner, we take Y = 6(2) /\/2, and enforce 
orthogonality to X just before normalization. 

It can be seen that the real and imaginary parts of 
(|3.8|) . used as X and Y, pass untouched through the 
Gram-Schmidt steps above. The choice we make of be- 
ginning with the original rfi instead will mean that our 
final TOqK differs from (|3.8|) by no more than a spin term. 
Unlike H3.8|l . the Gram-Schmidt procedure described can 
be used everywhere in the numerical domain, and guar- 
antees a smooth behavior moving between neighboring 
non-pathological and pathological points. 

Following the procedure as outlined here exactly re- 
produces the quasi-Kinnersley tetrad for Kerr-BL, and 
behaves well for Brill-Lindquist and Bowen- York-type bi- 
nary data. 



C. Spin-Boost Fixing and the Kinnersley Tetrad 

The tetrad obtained from this procedure will be trans- 
verse, and moreover, will be in the same transverse frame 
as the Kinnersley tetrad — differing only by a Type-Ill 
— or spin-boost - transformation. Such a transforma- 
tion will leave the scalar '^2 unchanged, but will have a 
strong effect on the radiative fields and ^'4, scaling 
and mixing polarizations. 

Lacking an unambiguous and natural way to lock down 
the spin-boost needed to obtain the Kinnersley tetrad 
from the quasi-Kinnersley frame member, we return to 
our knowledge of the Kerr background. It can be shown 
that for Kerr-BL, the null tetrad produced from the cor- 
rect transverse eigenvector (by following the construction 
in "41!) is 

fqK = ^i=[A,A,0,a] 



%/2AE 
1 



TOqK = 



V2AE 



[A,-A,0,a], 



(3.9) 



iasinO, 0, 1, — 



that is, the new tetrad is related to the Kinnersley tetrad 
fTT^ via: 



"■Kir 



^qK, 
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(3.10) 



This implies that the corresponding Kinnersley-tetrad ^'4 
can be obtained simply from: 



(3.11) 



Note that this transformation, unlike (|2.15|l for original 
Lazarus, does not involve any mixing of the Weyl scalars. 



D. Summary of the New Procedure 

Thus the new Lazarus procedure involves the following 
sequence of steps: 

(i) . construct at every point on the numerical grid a 

coordinate null tetrad of the form H2.9|l : 

(ii) . calculate the corresponding Weyl scalars the 

tetrad invariants X, J , and the speciality index 5, 
which we monitor during the full-numerical evolu- 
tion; 

(iii) . construct the matrix Q(i)(j) l|3.4|l . and find its 

eigenvalues and eigenvectors; 

(iv) . pick the eigenvalue A with the largest (complex) 

magnitude; 

(v) . use the corresponding unit-norm eigenvector V to 

form the principal (complex) spatial vector cr" = 

(vi) . construct a new null tetrad according (|3.t)|l . H3.7|l . 

and the pathology-avoiding steps following; 

(vii) . recalculate the Weyl scalars with this new tetrad; 

(viii) . determine the BL coordinates from the numerical 

ones via the transformations (|2.1(JI - [?.13|I : 

(ix). use the transformation 13.11|l to produce the 
Kinnersley-tetrad ^'4 and dt'^i] 



(x). evolve the Cauchy data -0 = / 
the Teukolsky equation (|2.4|) : 



^'4 and dti! using 



(xi). extract the gravitational radiation information, 
such as waveforms and total energy radiated at dif- 
ferent extraction times T. 



In the above, only steps (i)|(vii)|ca n be carried out in 
the full-numerical evolution; step |(ix)| is a post-processing 
step, as it still involves knowledge of the appropriate 
Kerr background parameters (M,a), and BL coordinates. 
While the new Lazarus approach outlined above does 
not solve this problem, it reduces the dependence of the 
Teukolsky data on the tetrad choices, by reducing the 
number of functions that depend on the BL coordinates. 



IV. NUMERICAL EVOLUTIONS 

We have applied the new transverse-frame-based tech- 
niques outlined above to three numerical regimes of in- 
creasing complexity: Bowen-York data for one spinning 
hole, Brill-Lindquist data for the head-on collision of two 
black holes, initially at rest, and transversely boosted bi- 
nary Bowen-York data, representing an "Innermost Sta- 
ble Circular Orbit" . Two of these data sets were already 
investigated in detail using the existing Lazarus approach 
in 13] ; our new techniques should allow some additional 
estimate on their validity. 

The results have been obtained from 3D evolutions per- 
formed with the standard ADM decomposition of Ein- 
stein's equations |60j |. together with zero shift and max- 
imal slicing lapse. It may be objected that this combi- 
nation of evolution systems and gauge conditions is not 
state-of-the-art. More sophisticated and stable systems 
are in wide use at present, and we would certainly expect 
them to produce much longer-lived stable full 3D evolu- 
tions. In particular, a new numerical-relativity frame- 
work (LazEv) has been successfully developed |i61i] and is 
currently available to the authors of this paper; it allows 
for evolutions of black-hole binaries using higher-order 
finite differencing, with the BSSN formulation of Ein- 
stein's equations and dynamical gauge conditions. Nev- 
ertheless, we have two reasons for using the "ADM -|- 
maximal slicing -I- zero shift" combination: First, this 
combination was that used for the original Lazarus tests; 
if we desire a fair comparison of old and new methods, fo- 
cused on the new tetrad methods we have developed, use 
of more modern and stable evolution methods will only 
obscure the results. Second, the introduction of more so- 
phisticated gauges will necessitate the reconsideration of 
the the numerical-to-Kerr-BL coordinate transformation. 
In particular, the use of other lapse evolution schemes 
may lead to a lapse shape significantly different to the 
maximal and Kerr-BL shapes. Also, the presence of a 
non-zero shift vector will alter the azimuthal-angle rela- 
tionship (|TT^ . 

We carried out the coding and testing of the concepts 
presented here using the Cactus framework. Post- 
processing of the 3D data and subsequent Teukolsky evo- 
lution was done with stand-alone codes. 

In these evolutions, we were also monitoring the ap- 
pearance of a single merged black hole using the apparent 
horizon finder AHFinderDirect .63.1 . 



A. Spinning Bowen-York Data 

The Bowen-York solution 0| representing a single 
black hole with angular momentum J should resemble 
— once the Hamiltonian constraint has been solved — a 
single Kerr black hole with gravitational radiation on top. 
The rotational symmetry of the solution about its spin 
axis (the coordinate polar or z axis) means that only radi- 
ation with TO = will be present, and that there should 



8 



be no net loss of angular momentum through gravita- 
tional radiation. Some analytic treatment of this data is 
presented in the Appendix Section 113 II based on Gleiser 
et al.'s perturbative analysis [65j| . 

In an early numerical investigation of the suitability 
of such data for Lazarus, we evaluated this data — us- 
ing a full-numerical solution to the constraints — and 
extracted the radiation content of the initial data via 
the new Lazarus method. We found that for small spins 
(J/M^ < 0.2), a spatial resolution of up to Af/80 would 
be necessary to identify unambiguously the leading-order 
real {£ = 2) and imaginary {£ = 3) parts of the Teukolsky 
function. Such a resolution is not feasible for a unigrid 
run, if the outer boundaries are to be outside the strong- 
field region. For this reason, we have used a larger initial 
spin for full-numerical simulations. 

This data was implemented on the initial numerical 
slice in Cactus, with bare mass m = 0.858 Madm, and 
spin angular momentum J = 0.553 M^^j^; after numer- 
ical solution of the Hamiltonian constraint, the ADM 
mass of the initial data was Madm — 1-0. To extend the 
physical domain, we used a "transition fisheye^ radial 
transformation of the numerical coordinates [llj. The 
explicit form we used in this paper can be found in Eqn 
(114) of 66J, with parameters {a = 7,s = 2,ro — 7); 
this produced a physical radial extent of 33.6 M from a 
numerical radial extent of 10.79M. We used three spa- 
tial resolutions — M/12, M/18 and M/24. Because the 
solution below the x-y plane is trivially related to that 
above the plane, we were able to evolve the upper bitant 
only. The 3+1 simulation died due to numerical insta- 
bilities before 25 M of coordinate time at the coarsest 
resolution, before 19 M at the medium resolution, and 
before 18 M at the finest resolution. 

For such a weak-field case, we expect the Lazarus tech- 
niques to be applicable from a very early time. An ap- 
parent horizon is already present in the initial data, and 
Fig. n shows the value of the (real part of the) speciality 
index S along the x axis outside the horizon, evaluated 
at several times in the evolution. At all times after the 
initial slice, the deviation from the Kerr value of 5 = 1 
is at most ~ 2%. Thus the spinning BY data should 
certainly be within the linearized regime almost imme- 
diately. (Nevertheless, Lazarus results will be sensitive 
to how well the numerical space-time satisfies the coor- 
dinate assumptions of Section FlI B\ ) 

We analyzed the evolved data with both the original 
and new Lazarus approaches. At every M of coordi- 
nate time evolution, the Weyl scalars were calculated 
using both the original and new tetrads, then mode- 
decomposed and saved to file. We post-processed this 
data to produce the Teukolsky Cauchy data (i/;,9ti/'), 
which was then evolved using a code that solves the 
Teukolsky Equation. The total extracted energy as a 
function of extraction time T is shown in Fig. |3 Be- 
tween extraction times of IM and 5M, we see an emitted- 
energy plateau at around 1.25 x 10~^ M, common to both 
original and new Lazarus, with a deviation of ~ 3.5% be- 



Spinning BY hole, Re(S) along x axis 
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FIG. 1: The real part of the S invariant along the x axis for 
Spinning BY data, at different extraction times T. 
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FIG. 2: Total radiated energy for Spinning BY data as calcu- 
lated from original (LI) and new (L2) Lazarus methods, as a 
function of extraction time T. 



tween old and new values. 

We would expect that once we have entered an energy 
plateau, we should remain there, at least until the radi- 
ation has left the outer boundary of the full-numerical 
domain. However, for later extraction times, from 6M 
to lOAf , we see a monotonic decrease in emitted energy, 
leading to a low point of ~ 0.88 x 10~^M, a drop of 
^ 30 %. This degradation is common to both original 
and new Lazarus. 

This behavior may be an artifact of the underlying 
ADM evolution system, or other aspects of the Lazarus 
procedure; even for the relatively high spin chosen here, 
spinning Bowen-York is a weak source of gravitational 
radiation, and presumably extremely sensitive to details 
in the simulation. Further investigation with newer, more 
stable and accurate, evolution systems may illuminate 
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Spinning BY hoie, m=0, Lazarus 2 
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FIG. 3: Evolution of the m = Spinning BY waveform for 
the new (L2) Lazarus method, evaluated at r« — 30 M along 
maximal amplitude directions: 9 = 90° for J?e(*I'4), = 45° 
for Im{^4). The five curves represent five different extraction 
times in the energy plateau of Fig. |21 



the situation. 

The level of emitted energy is consistent with the full- 
numerical results shown in |67j and the perturbative 
results of Gleiser et al. : they obtain a total emitted 
energy of ^ 2.5 x 10"'' M for a spin of Jadm = 0.55 A-P. 
Although significantly larger than our peak energy of ~ 
1.25 X 10^* Af, it falls into the "factor of two" accuracy 
the authors assess for their method at high spins. 

Selected Teukolsky-evolved waveforms (the rescaled 
Weyl scalar ^'4), measured at (r* — 30 M, 9 = 0), are 
shown in Fig. |31 The waveforms shown are for the new 
tetrad method only (old-method waveforms are indistin- 
guishable). There is good agreement especially between 
extraction times T = 3M and 5M. For later extraction 
times, the waveforms lose coherence after a few wave- 
lengths. 




T" 



— T = OM 
■■■■ T = 2M 
-- T = 4M 
■-■ T = 6M 
■-■ T = 8M 
-■■ T=10M 

— T=12M 



z/M 



FIG. 4: The S invariant along the z axis for Head-On data, 
at different extraction times T. 



B. Brill-Lindquist Data 

This test revisits the case of head-on collisions of Brill- 
Lindquist data, treated previously in the Lazarus context 
[13, Data of this type was studied in a perturba- 
tive form in Ref. [H3 i based on the full- numerical results 
available at the time, perturbative results overestimate 
the total emitted energy for large separations. Quali- 
tatively similar conclusions can be drawn using Misner 
initial data Q. 

The specific initial data chosen has two holes with 
equal bare masses m = 0.50 M, with a coordinate sep- 
aration along the y axis of 2.303 M, yielding a proper 
horizon-to-horizon separation of L = 4.9^, and no lin- 
ear or angular momentum (this data set was referred to 
as "P = 0" in 13J). The symmetry of the problem al- 
lows us to use only the first octant of the numerical grid 
during evolution. Otherwise, the grid extent and fisheye 
transformation were identical to those for the spinning 
Bowen-York case above. 

Again, we evolved the initial data using the "standard 
ADM" evolution scheme with maximal slicing and zero 
shift. The finest resolution run crashed due to numerical 
instabilities before 20Af of coordinate time. Neverthe- 
less, a common apparent horizon was found at coordi- 
nate time t « 8M, and there is reason to believe that a 
common event horizon is present from t fa 3 M |69l | and a 
continuous potential barrier surrounding the strong-field 
region around even earlier. We show the evolution of the 
speciality index S in Fig. 01 While deviations from unity 
are always large near the punctures, the overall devia- 
tion outside the eventual common horizon location have 
dropped to below 10% by T w 8 M, indicating we have 
entered the linear regime. 

Fig. [Slshows the total radiated energy from the Teukol- 
sky evolution, as a function of the extraction time T for 
original and new methods, for the m = (upper panel) 
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FIG. 5: Radiated energy for Head-On collision as calculated 
from original (LI) and new (L2) Lazarus methods, as a func- 
tion of extraction time T. Note the extension of the plateau 
with increasing resolution. 



TABLE L Comparison of emitted E for Head-On collision. 
Direct results {r,^obs = 7.46 M) are from 61]. E and ctb are 
the mean and standard deviation from four plateau values. 



Source 


\E\/M 


as/M 
xlO-'' 




Direct 


6.6 




LI (8 


- IIM) 


6.827 


0.446 


L2 (8 


- IIM) 


6.811 


0.179 



Head-On collision, m=+2, Lazarus 2 
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FIG. 6: Evolution of the m = +2 Head-On waveform for the 
new (L2) Lazarus method, evaluated at (r, = 3QM, 6^ = 0). 



ties of the method when nonlinearities are present the 
agreement at later times shows when a safe, linear regime 
is reached. 



and m = +2 (lower panel) modes. For both sets of 
modes, original and new Lazarus deviate at early extrac- 
tion times. While the original Lazarus results reproduce 
what was seen in Fig. 5 of [l^l — the radiated energy 
reaches a crude plateau between extraction times of 4M 
and lOM — the new-tetrad-produced energy reaches a 
level plateau only after 7 M. After this time, both meth- 
ods agree well until 10 M, but with a noticeably more 
level plateau in the weaker m — curve for the new 
tetrad. 

Table m compares the total energy calculated with old 
and new Lazarus with a direct-extraction energy taken 
from 'oil (this extraction was from a rather close detector 
position, r^^^obs ~ 7.46 M). While both Lazarus figures 
differ from the direct result, the new Lazarus plateau - 
measured from T — 8M and T = 11 M ~ has less than 
half the standard deviation of old Lazarus, indicating a 
more stable plateau. 

Fig. shows the dominant-mode (m = -1-2) Lazarus 
waveform for the Head-On collision, measured at (r* = 
30 M, 9 = 0). While the difference between curves for 
early extraction times can give an idea of the uncertain- 



C. QCO Data 

The final test of our method again revisits initial 
data treated in ■ Bowen-York binary black- hole data, 
where the holes have zero spin, but are boosted in a di- 
rection transverse to their separation, to achieve a net 
orbital angular momentum J. 

Khanna et al. [t^ investigated such data to first per- 
turbative order, through Zerilli and Teukolsky schemes. 
The emitted energy curves from the two methods diverge 
at the level J /NP « 0.4—0.5, and this the authors take as 
the limit of applicability of linear perturbative theory. At 
the highest analyzed spin, J/M^ = 0.6, the Teukolsky- 
and Zerilli-calculated energies are ~ 8 x 10^"^ M and 
~ 5 X 10"'^M respectively. 

Here, we use the binary Bowen-York data with equal 
bare masses m — 0.45 M, located at coordinate positions 
y = ±1.1515 M, and with Bowen-York "boosts" P = 
±0.335 M in the x direction. These result in a physical 
throat-to-throat separation L = 4.9 M, and a total ADM 
angular momentum J = 0.11 NP. We refer to this data 
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QCO, Re(S) along x and z axes 
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FIG. 7: The S invariant along the x (upper) and z (lower) 
axes for QCO data, at different extraction times T. 

set as "QCO", following the classification scheme of [T3| : 
it was originally suggested by Baumgarte [tJ (adapting 
Cook's "effective potential" method "7^ to punctures) as 
a model for the so-called Innermost Stable Circular Orbit 
(ISCO) of two equal-mass black holes. 

In contrast to the two previous test cases, numerical 
instabilities will kill the QCO data evolution before a sin- 
gle common apparent horizon has appeared (though a 
common event horizon may already be present). Thus 
Lazarus extraction and evolution of Tcukolsky data 
would not appear to be completely justified at any time. 
Nevertheless, [13| found a small plateau in the emitted 
energy, indicating a range of extraction times where the 
system has effectively linearized prior to the simulation's 
crash. 

For this data, the Lazarus procedure was carried out 
initially using a zeroth-order parameter estimate taken 
from the ADM energy and angular momentum of the 
initial data: M — 1, a — 0.77. However, since the emit- 
ted energy and angular momentum were at the level of a 
few percent of the total, the subsequent drop in back- 
ground mass and spin might be significant. For this 
reason, we iterated the method with suitably reduced 
background mass and spin parameters: M = 0.974, and 
a = 0.675/M = 0.693. 

To demonstrate how quickly the system appears to re- 
lax to Kerr, we show in Fig. [7|the S invariant along the 
X and z axes. Along the z axis, S clearly takes time to 
settle down, only uniformly falling within [0.5, 2.0] rather 
late in the simulation, after T = 8 M. 

Fig. |S1 shows the emitted energy and angular mo- 
mentum, respectively, for both original and new Lazarus 
methods (the dominant m — +2 mode waveform is plot- 
ted for old and new Lazarus in Fig. and Fig. 
respectively). The original Lazarus results can be found 
in Fig. 22 of 0| (note that the latter figure contains the 
total energy summed across m modes). There is signifi- 



TABLE II: Comparison of emitted E and J for QCO merger. 
Direct results (rt^obs = 14.03 Ai") are from 73j. E, J are the 
means and oe, crj the standard deviations from five plateau 
points (indicated in first column). 



Source 


\E\/M 
xlO"^ 


<te/M 
xlO"^ 


|J|/A/2 
xlO"^ 


xlO"^ 




Direct 


2.8 ±0.2 




11.0 ± 1.0 




LI (9.0 - 


11. OA/) 


2.570 


0.067 


9.493 


0.434 


L2 (8.5 - 


10.5A.f) 


2.422 


0.031 


8.553 


0.175 



cant deviation between original and new Lazarus results 
up to the extraction time of T « 6 M. The two methods 
agree well for T € [7 M, 8 M] , but begin to diverge again 
around the time new Lazarus develops plateau values of 
around 1.1 x 10"^^/ ^nd 4.25 x 10"^ j\,/2 f^j. emitted 
energy and angular momentum, in the mode m = -1-2, re- 
spectively. Working with the more time-resolved data of 
the inserts in Fig. O we have identified the plateaus for 
old and new Lazarus, and calculated the means and stan- 
dard deviations. We present these in Table UTI along with 
direct-extraction results (r^^obs = 14.03 M) from [73|| . 
While both old and new means differ from the direct- 
extraction figures, new Lazarus's plateau is flatter. 

Because the QCO data takes time to plunge to a single 
hole, the onset of linearization will be delayed relative 
to, for instance, the head-on case. The fact that the ob- 
served plateau does not persist for longer than 2 M is 
a consequence of the instability of the the simple combi- 
nation of "ADM -f maximal lapse -I- zero shift" we have 
used in our full-numerical evolutions here. In this sense, 
QCO -I- ADM is a "marginal" Lazarus case; to establish 
the plateau unambiguously will require evolutions that 
last for several more M. This will necessitate more ad- 
vanced evolution systems and gauges, and is addressed 
in the Discussion below. 



V. DISCUSSION 

The original Lazarus method was a successful synthesis 
of full-numerical and perturbative methods in numerical 
relativity. It was responsible for the production of the 
first waveform for black-hole binary mergers and close 
orbits, before these regimes were accessible to a complete 
full- numerical treatment. 

We have presented in this paper the first in a series 
of proposed improvements to the Lazarus method — im- 
provements which will, we hope, make Lazarus more am- 
bitious in its reach, and more rigorous in its grasp. The 
developments herein have focused on the important area 
of tetrad determination. We can now obtain — in the 
"real time" of the full-numerical evolution — a trans- 
verse tetrad that will differ from the Kinnersley tetrad 
only by a spin-boost. 

We described in Section IIIII how this new method 
of tetrad determination can be incorporated into the 
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FIG. 8: m — +2 radiated energy and angular momentum 
for QCO data as calculated from original (LI) and new (L2) 
Lazarus methods, as a function of extraction time T. 



FIG. 9: The two polarizations of the m = +2 QCO wave- 
form for the original (LI) Lazarus method, evaluated at 
(r. = 31M,e = 0). 



Lazarus project, eliminating the artificial mixing of 
monopolar, longitudinal and transverse modes that was 
necessary in original Lazarus. To test this new aspect, we 
have applied it on the same footing as original Lazarus 
— to the short-lived ADM evolutions of three different 
types of black-hole data. 

Our numerical results with this method, reported in 
IIVI have allowed us to test this mode mixing for the first 
time. The good agreement between old and new methods 
for the single Spinning Bowen-York and Head-On cases, 
where linearization is seen well before code break-down, 
can be seen as a validation of the applicability of the orig- 
inal Lazarus mixing formulas; additional evolution time 
may be needed to resolve the marginal case of QCO data 
satisfactorily. In contrast, the areas of difference — seen 
for instance, before T = 7 M in Fig. Ofor the Head-On 
case — indicate better where our assumption of linear de- 
viation from Kerr may not have been justified. With this 
alternative path to the Kinnersley tetrad in Lazarus, we 
have produced error estimates for the tetrad construc- 
tion procedure; we can estimate the contribution that 
errors in this procedure make to the overall uncertainty 
in Lazarus energies. It should be noted, however, that 



this error is unlikely to dominate the total uncertainty 
in emitted energies. Comparison of Lazarus (both old 
and new) results with more direct extraction techniques, 
shown in Tables ^ and ^ indicate that other important 
factors must still be addressed. A reasonable estimate 
of theoretical error in waveforms will be needed for their 
application to analysis of gravitational-wave detector sig- 
nals in the future. 

For the next stage in the Lazarus project, we plan to 
use modern evolution schemes, such the one provided by 
the LazEv framework j^, which are more stable and 
accurate than ADM. LazEv currently supports higher- 
order finite differencing for the BSSN formulation of Ein- 
stein's equations with the choice of several dynamical 
gauge conditions, and should allow for extended energy 
plateaus, extending from the time of first linearization of 
the system until the time the radiation leaves the system. 
This added full-numerical evolution time will aid in the 
unambiguous identification of the linearization time for 
marginal cases such as "QCO" above. 

At present, it is unclear how to relate Boyer-Lindquist 
time to the numerical time developed with standard nu- 
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FIG. 10: The two polarizations of the m — +2 QCO wave- 
form for the new (L2) Lazarus method, evaluated at (r, = 
31 M, 8 = 0). 



merical gauges. Maximal slicing was discussed in Section 
III El future refinements of this may use analytic insi ghts 
into the late-time shape of the maximal lapse (see [tJ 
for work on maximal slicings of Schwarzschild). Com- 
monly used dynamic slicings (e.g., "l-|-log" slicing) en- 
sure a lapse function with the same qualitative shape, 
falling smoothly off to unity at large distances. The exact 
shape will vary greatly in the near-field region, however, 
and this will affect both old and new Lazarus; this is an 
interesting problem, which we hope to address in future 
work. 

Treatment of the shift is decoupled from the lapse in 
the Lazarus approach; the shift correction via adjust- 
ment of is independent of the slicing analysis. Co- 
rotating coordinates can already be accommodated with 
this method. We expect to treat more complicated evo- 
lution shifts in a similar way. 

The quasi-Kinnersley frame has further useful appli- 
cations. For direct radiation extraction at a finite ob- 
server location in the 3D numerical grid. It also pro- 
vides a closer interpretation in terms of radiation for 3- 
dimensional visualizations. Originally one uses the nu- 



merical tetrad to evaluate ^'4. In order to give this a 
direct interpretation in terms of radiation, one implic- 
itly assumes that the observer is sufficiently far from the 
strong-field region that the background space is almost 
flat, and the numerical tetrad is a good approximation 
to Kinnersley. 

With the quasi-Kinnersley frame, we can get much 
closer to the true Kinnersley form at finite observer lo- 
cations. This should supply us with a more robust wave- 
form at distant locations, and allow us to approach the 
near-field zone with greater confidence. While imprac- 
ticable for short-lived ADM evolutions, modern evolu- 
tion systems should produce several periods of directly 
extracted waveforms at different extraction radii. Outer- 
radius waveforms may be compared with Lazarus results; 
inner waveforms may be compared with outer ones as a 
test of the quality of strong-field waveforms in general. 

Aside from the longer full-numerical evolution times, 
we plan to investigate further several issues to improve 
the generality of Lazarus. 

The transformations used to extract Boyer-Lindquist 
coordinates from the numerical ones perform well, but 
there is much room for improvement. One obvious — 
and theoretically well-founded — step that can be taken 
is relaxing the assumption that there is a one-to-one re- 
lationship between the numerical radial coordinate i?nuin 
and the Boyer-Lindquist r. Since we know the theoreti- 
cal form of the invariant X as a function of (r, 9) , we can 
imagine using X everywhere to derive both r and 9. We 
have performed initial investigations of this possibility, 
but found that the off-equator radial transformation de- 
rived does not reach sufficiently close to the horizon, at 
least for the most difficult cases, such as the QCO. Longer 
full-numerical evolutions may improve the performance 
of this technique, as a horizon appears and circularizes. 
We also plan to revisit the choice of the time slicing as 
this might be crucial to process longer-term fully nonlin- 
ear evolutions. 
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APPENDIX A: NUMERICAL-KINNERSLEY 
TRANSFORMATIONS FOR KERR DATA 

In Boyer-Lindquist (BL) coordinates, the numerical 
tetrad — defined by H2.9|) with orthonormalized spher- 
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ical coordinate directions for the e 



(i) 



takes the form: 



AE' 



'A 2aMr 



n A 2aMr 



0,0, 



1 



(Al) 



^Kin 



[i^(D- 1)* 



num 





-r2(i?- 1-2^2)*™ 

2(Z? + l + 2v42)<lf™™ (A7) 



It can be easily verified that these transformations, ap- 
phed to the numerical-tetrad Weyl scalars of Kerr given 
in ljA2|l . produce the Kinnersley values: 



^f^2 = 0. 



Such a tetrad will differ strongly from the Kinnersley 
tetrad; as a consequence, all Weyl scalars calculated from 
it will be non-zero. For Kerr-BL, these values will be: 



*2 



M 
M 
M 

2nc^ 



[3iAVA2^]. (A2) 



^(3A2-f7)]. 



For Kerr-BL coordinates, the numerical null tetrad 
(|A1I) used in 3-1-1 numerical calculations can be trans- 
formed to the Kinnersley tetrad 1)2. 7|l via a combination 
of null rotations and spin-boosts: 



h 



"Kin = 



Fa 
2 



{in 



1)1. 



p-1 



2 



0}, 

+ {D 
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-iA{l. 



Note that this is a corrected form of the transformations 
that appeared in Eqns (5.9a-c) of Paper I. 

The dependence of the Kinnersley-tetrad on the the 
numerical-tetrad scalars was given in (|2.15() . for pertur- 
bations of Kerr-BL coordinates. The following are the 
corresponding expressions for the remaining Kinnersley 
Weyl scalars (the mixing function A and spin-boost func- 
tions Fa and Fb are given by H2.14|l . and D = y/A^ + 1): 



Kin 



(A4) 



Kin 



Kin 



[{D + lf + 4iA{D + l) 

-6 A^ -4iA{D-l) ^-^"^ 
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[-i A{D + 1) ^g"'" + 2{D + 1 + 2A^) ^-r'" 
+6iAD + 2 (D - 1 - 2 yl2) (A5) 
^iA{D-l)^r"']FAFB/A, 

4-5""^ -AiAD ^-^""^ + {A + 6A^) *7"^-H 
" /4, (A6) 
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APPENDIX B: PERTURBATIVE RESULTS 
1. Weyl Scalars for Spinning Bowen-York Data 

Gleiser et al. investigated the low-spin behavior 
of the Bowen-York data, obtaining an analytic solution 
for the conformal factor accurate up to O(e^), where e = 
J/M^, and evolving the extracted radiative modes via 
the Zerilli formalism. 

For zero spin, the Bowen-York solution reduces to the 
Schwarzschild solution in isotropic coordinates, and is 
automatically constraint-satisfying; for small spins, the 
dominant radiation {£ — 2) should scale as e^, while the 
next mode {£ = 3) will scale as e"^. 

As an estimate of the radiation content of this data, 
we may calculate the Weyl scalars of the approximate 
solution above; using the numerical tetrad (|Aip . we find 
on the initial slice: 



*2 



1536 A/3 i?5 ^ ^ . „ , . 2 ^ 

5 {2R + My^^ ' 
_ 384 ]VP e 
-fs - ~(2i? + M)i5 

64 M 



{iAi+eBi+ie^ Ci), 



5 (2i? + M)i4 
+1 ^2), 



(A2 + i e S2 + C2 



(Bl) 



where the e coefficients in parentheses are: 

QAMR + M^) 



Bo 
Al 
Bl 
Ci 

A2 
B2 
C2 



D2 



M 



~{2R + Mfi-iR^ 
-96 M^i?^ cos 61, 
-5(2i?- Af) (2i?-f Af)^ 
16 M R{2R + Mf {2R- M) cos 9, 
-IQM^{2R + Mf (14 i?2 - MR + 
+672 Af3 R^{2R- M) s\t? 9, 
5{2R + Mf, 
mM R{2R + Mf cos 9, 
12 M^^ {2R + Mf 

X [-{2R + Mf {8R^ - iM R~ M^) 
+A R^ (12 i?2 - 4 Af i? + 3 Af2) sin^ 9] 
192 Af* i? cos 6* 

X [{2R + Mf (4 i? + Af ) - 36 R^ sin^ 
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For simplicity, we factor out k = 64 M R^/{2R + Mf 
from the Q matrix; to 0(e), the three eigenvalues of the 
reduced Q matrix here are: 



Ml 



2 + ie 

-1 - ? 
-1. 



24 M R cos e 
(2i? + M)2 ' 
24 Af i? cos 61 
' (2i? + Af)2 ' 



The desired eigenvalue is /ii, which tends to —2 x /i2 or 
/i3 as e — > 0. When the factor n is put back in, the full 
eigenvalue is 



Ai = 128 



M R%2R + Mf + UieMR cosO] 
(2 i? + M)8 
Af (r + 3 i a cos 6*) 



where we've interpreted a = A/e as the (dimensional) 
Kerr spin parameter, reintroduced the Kerr-BL radial 
variable r = R + M + {AP - a^)/{AR), and neglected 
terms of higher than linear order. The corresponding 
eigenvector will be (any multiple of): 



V 



4iM Rsm9{2R - M) 
(2i? + Af)3 



e,0,l 



Note that V becomes real at 6* = 0, and a.t R — A//2. 



2. Weyl Scalars for Head-On Data 

The data used for the head-on run was of the Brill- 
Lindquist type, with holes separated along the y axis. 

As an estimate of the radiation content of this data, 
we may calculate the Weyl scalars of the approximate 
solution above; using the numerical tetrad HA1|I . we find 
on the initial slice: 



-48 M 



+96 i e-" M' 



48 Af^ 



3 (COS^ (f> — COS^ 



(2i? + A/)5 
cos 6 sin (j> cos (p 
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R + M) sin 6 cos 6 sin^ . 
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+4:8 ie^ ■ 
64 M i?3 
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\ R + M) sin sin cos < 



-I- 16 e-" 



{2R + Mf 
Af3 



{2R + Mf {2R + My 

x(24 R^ + 2M R + M'^)T{e,(j)), 



(B2) 



where we define T{6, 0) = 2 — 3 (cos^ 9 + siif 9 cos^ 4>). 
For simphcity, we factor out k = 64 A/ R^/{2R + Mf 
from the Q matrix; the three eigenvalues of the reduced 



Q matrix are: 

Ml = -1, 
Ai2 = 2 -he 

Ma = -1 - e 



A/2 (24 R^ + 2MR + M'^) T{9, (f) 

2R^{2R + M) 
2 M^ {24: R^ + 2 M R + M^)T{9, 4 
2R^{2R + M) 



Here the principal eigenvalue is obviously /i2, and the 
rescaled equivalent is: 



128 Afi?^ 
{2R + Mf 

2 32 Af3 R^{24R^ + 2MR + M^) T{9, < 



i?3 {2R + My 



128 MR^ 



{2R + Mf 

32 (P M R^ (24 R^ + 2MR + M^) T{9, </>) 
^ i?3 {2R + My 

with a related normalized (to 0(e^)) eigenvector 



V = 



-e2 Af 2 (8R + Af ) sin 6* cos 9 sin^ 0/2 f?^ 
e2 Af 2 (8 f? -I- Af ) sin 9 sin cos (/)/2 f?^ 
1 

-(P{8R + Af ) sin 6* cos 9 siif <j)/2R^' 
(f{8R + M) sin 9 sin (/) cos (/)/2 f?^ 
1 



Note that in this case, the normalized eigenvector is 
manifestly real everywhere, and so would lead to a de- 
generacy problem in the reconstruction of the quasi- 
Kinnersley tetrad, of the type described in Section IlIIBI 



3. Weyl Scalars for QCO Data 

The data used for the QCO run was of the binary 
Bowen-York type. We can try to determine some prop- 
erties of the slow-close limiting form of this data. This 
was addressed by [t^I, who treated the binary data — 
with zero-spin holes separated by L in coordinate space, 
and transversely boosted with momenta ±P — as a per- 
turbation of Schwarzschild with perturbation parameter 
e = LP/M^. The authors worked solely in the Zer- 
illi formalism (they point out that such data is only a 
perturbation of Schwarzschild, not of Kerr), and only to 
0(e), where the Hamiltonian constraint did not need to 
be solved for consistency. 

Taking this data (separation along the y axis, and 
boost in the x direction), the numerical tetrad yields 



'I'2 

^-1 



64 Af R 
(2f? + Mf 

384 e Af 2 f?4 sin 9 



[{2R + M f - UieMR cos 9] 



(2f? -I- Af )9 



[4(2f? + Af ) cos sin (j) cos f 
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-6iR-iM + 4:i{2R + M)cos^4>] , 
IQ^eAP R'^(2R- M) 
" {2R + Mf 

X [cos0 sin(^(l + cos^ 9) — icos6{l — 2cos^ ^)] . 

For simplicity, we factor out k = Q4MR^/{2R + Mf 
from the Q matrix; the three eigenvalues of the reduced 
Q matrix are: 



lJb\ = 2 — ie 
Hi = -1 + ie 



24 MR cos e 



{2R + MY ' 
24 M R cos e 
{2R + Mf ' 



/"3 



-1. 



The desired eigenvalue is /^i, which tends to — 2 x ^2 or 
/X3 as e ^ 0. When the factor k is put back in, the full 
eigenvalue is 



Ai 



128 



MR%2R + Mf -UieMR cos 6*] 
{2R + Mf 



M{r + 2>ia cos 6) 



where we've interpreted a = —Me as the (dimensional) 
Kerr spin parameter, reintroduced the Kcrr-BL radial 
variable r = R + M + (M^ - a'^)/{4:R), and neglected 
terms of higher than linear order. The corresponding 
eigenvector will be (any multiple of): 



V = 



AiM R sine {-6 R - M + 4{2R + M) cos^ cj)) 

(2i? + Mf 
16 z M i? sin cos sin cos 0) 



{2R+Mf 



Note that V becomes real for certain angular posi- 
tions: e.g., ^ = for all (j), and ^ = 7r/2 for = 

isJ^m)^^/^^ovR»M. 
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